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The Feynman Calculus 


This is chapter 6 in Griffiths. 
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The Feynman Calculus 


We are now going to skip (for now, will cover it later) chapter 5 on 
"Bound States" and go straight to chapter 6 of Griffiths. 


Chapter 6 introduces the Feynman Calculus with an excellent toy 
theory before we apply the methods to QED, QCD or electroweak 
theory. 


Jyothsna (CHEP) Nuclear & Particle Physics Course October 2018 Lectures 309 / 431 


Decays and Scattering 


@ We have three experimental probes of elementary particle 
interactions: 
> bound states (covered later) 
> decays 
> scattering 
9 Nonrelativistic quantum mechanics (in Schrodinger's formulation) 
is particularly well adapted to handle bound states 


© By contrast, the relativistic theory (in Feynman's formulation) is 
especially well suited to describe decays and scattering. 
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Observables 


We want to relate experimental measurements to theoretical 
predictions. 


e Decay widths and lifetimes 
T = h/r (units of energy) 


e Scattering cross-sections 


G is the total cross section 


do is the angular distribution 
v is the energy distribution 
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Decay Rates 


The decay rate F is the probability per unit time that any given particle 
will disintegrate. So, for a group of particles we have 


dN = -T Nat 
N(t) = N(0)e"'' 
The mean lifetime is defined as the reciprocal of the decay rate 


=i 
fer 


Most particles can decay via multiple distinct channels and the total 
decay rate is the sum of the rate in each channel 


n 
[tot = X li 
ici 
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Breit-Wigner Resonance 


Wavefunction for a particle with rest energy Ep and width T is 
V (t) - V(0)e- En —-ir/2)t 


Fourier transform x(E) x TER PJA 


Breit-Wigner formula : Ix(E)|? « UR EUGTUA 


The production cross section (rate of production per incoming 
particle) is described by the Breit Wigner resonance formula 


r? 


o (E)“ ——— 
(E - Mc + 02/4 


where M is the central 
mass of the particle and 
T is its width, 
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Decay Rates 


The Branching Ratio is the fraction of all decaying particles which 
decay via a specific channel and is logically defined as 


li 


Bl 
O Tiot 


OK, those are the important terms for decays, how about scattering? 
Well, ultimately we are interested in the cross section...but need to 
develop some formalism to get there... 
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Scattering 





A cross section is a geometric idea you can relate to from macroscopic 
collisions. We need to generalize it. 
Q The interaction between the projectile and the target can be 
"long-range". 
Q The cross section is not a sole property of the target but rather a 
joint characteristic of the projectile and the target. 
Q There are inelastic processes in which the final state particles are 
different from the initial state particles. 
Just using classical mechanics we can describe scattering as a way of 
relating an impact parameter to a scattering angle. 








Scattering center 
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Hard-Sphere Scattering 


Example 6.1 
A particle bounces elastically off of a sphere of radius R. 








Figure 6.2 Hard-sphere scattering. 
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Hard-Sphere Scattering 


b = Rsina, 
2a +0 = m 
; ft 8 0 
sina = sin(5 — 5) = cos(5) 
=> b = Reos(5) 


This is the relation between 9 and b for classical hard-sphere 
scattering. 
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The Differential Cross Section 









8 
A 
NN 


= 









If a particle passes through an infinitessimal area do it will be deflected 
into a corresponding solid angle dO. The larger we make do the larger 
is dQ and the differential cross section is 


do = |bdbdd¢|,dQ = |sin 0 dé do| 
sa = lina (ao) 
dO — {sin@ \da 
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Short-range interaction: Hard-Sphere Scattering 


Examples 6.2 and 6.3 


Find the differential and total cross sections for hard sphere scattering 
with a sphere of radius R. 





From example 6.1 we know that b = R cos(0/2) so 
db R in (2) 
da 242 


do  Rbsin(0/2)  R?cos(0/2)sin(0/2) | R“ 
dO . 2snO 2 sin 0 4 
Giving a total cross section of 


2 
o= [ gan - [ an =n 
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and 





Long-range interaction: Rutherford Scattering 


Example 6.4 - Rutherford Scattering 
A particle of charge qı scatters off a stationary particle of charge qo. | 
In classical mechanics the formula relating b to 0 is 


„MR 
b = SE cot(0/2) 


So, the differential cross section is 


(= (z) 

sin 0 do 4E sin?(0/2) 

The total cross section is actually infinite... 

This is related to the fact that the Coulomb potential has infinite range.. 


do 
dQ 
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Fermi’s Golden Rule 


Fermi’s Golden Rule 


2 
transition rate = X WP x (phase space) 





e The amplitude M contains all of the dynamical information. Use 
Feynman diagram/rules to calculate this. 

e The phase space is a kinematic factor. The bigger the phase 
space the larger the transition rate. 

© Alternate terminology: 


> Amplitude « Matrix Element 
> Phase Space o Density of Final States 


Jyothsna (CHEP) Nuclear & Particle Physics Course October 2018 Lectures 321 / 431 


Fermi’s Golden Rule for Decays 


Golden Rule for Decays 
Suppose particle 1 decays as 


152+3+4++n 


The decay rate is given by 
r= ges. | Pleny“ (p — p — py — Pn) 
^ 2hm E c Pn 


x [ons (rj — mid) 0 (W) 
j= 


d'pj 
(27)“ 











Where the ó is a Dirac-Delta, 8 is step function, p; is the 4-momentum 
of the i-th particle, the decaying particle is at rest (mc, 0) and S is 1/j! 
for each group of j identical particles in the final state. 


Jyothsna (CHEP) Nuclear & Particle Physics Course October 2018 Lectures 322 / 431 


Dirac Delta & 0 functions 


Dirac delta function 


: oo 
i-a-|2 ias and [ sea 


1 5(x-a) 
| 


f(x)š(x — a) = f (a)à(x — a) 


[fen — de fia 








Step function 


CIN (x « 0) 


1 (x>0) 
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Fermi’s Golden Rule for Decays 


The dynamics of the process is contained in amplitude, M, which we 
will calculate later using appropriate Feynman diagrams applying 
Feynman rules. 

The rest is phase space which tells us to integrate over all outgoing 
four-momenta, subject to three kinematical constraints: 


1. Each outgoing particle lies on its mass shell: p? = mic? 
(which is to say, E? - Pj c? = mic*). This is enforced by the 


delta function à (7 — mid), which is zero unless its 


argument vanishes. 

2. Each outgoing energy is positive: pj — Ej/c » 0. Hence the 8 
function. 

3. Energy and momentum must be conserved: py = p2 + pi 
+ pn. This is ensured by the factor 8“ (py — p2 — ps +++ —pn). 
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Fermi’s Golden Rule for Decays 


e Thumb rule for factors of 27x: 
Every 6 gets (27) and d gets 1/27. 

e Four-dimensional 'volume' elements can be split into spatial and 
temporal parts: dtp = dp?d?p 

9e Using the properties of ô function: 


5 (p = mc) =6 [p° = p = mc] 


8(x* — a’) = x [5(x — a) + 5(x+a)] (a> 0) 


1 
0 042 2 2,2] _ o0 | 2 2.2 
9 (p°) s [(p") -p ome) = z = zr yP m) 


0 function kills the spike at p? = — yp? + m?c? 
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Fermi’s Golden Rule for Decays 


Thus the total decay rate is given by: 


S 
— ^ 


o 1 d dp 
ne s +m 22 2 [pet ma (2x)? 


with 


Pp yp + mpc? 
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Two-particle Decays 


With only 2 particles in the final state we have: 


S 8*(p — p — p3) 3 3 
T ===- 2 
32n thm, [^ d'p,x d'p, 


The four-dimensional delta function is a product of temporal and 
spatial parts: 


8*(p — p2 — ps) = 8 (pi — p? — p3) 5° (p — P2 — Ps) 
With particle 1 at rest, the decay rate is: 


— ien cuin piede) 
— 39-2komw. — 
327x2hm [p2 + mic |p + mic 


x 8? (pa + ps) d?p, d?p, 
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Two-particle Decays 


ps integral is now trivial: in view of the final delta function replace 
ps > -pe 


ur me vte- fee) 
VPs + mic? pi + mie ý 


Using spherical co-ordinates p» — (r, 0, p) and d?p» — r“sinOdrd9dý 
with r = |pal: 


^ 32z?hm, 


s (me- fremde - fremd 


S 2 
— 32x?hm, J M 
x 1? sin drd9 dý 
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Two-particle Decays 


M is only function of r and integrating over the angles 
J sin @ d0 dý = 4n gives 


5 (me- / + mie - e + mie) 
Jr + mier + mic? 


r= ze | |.4(r)| r dr 


Let 
u = fret mic? + fr? + mic? 
so 
du — ur 
dr JP + mig? + mic? 
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Two-particle Decays 


rx] lA (r)|28 (mc — u) < du 
— 8rhm, J, (r) : u 


m» --ma3)c 


Last integral sends u to m,c and hence r to: 


c 
fg = zm V" + m3 + mí — 2mm — 2mim3 — 2mm 
1 
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Two-particle Decays 


myc = \/r2 + m2c? + \/r2 + m2c2. Square: 
1 2 3 q 








mic =r +m? +r 2 + mac? + 24/72 + m ejr +m c 





c . 
3 (mj — m$ — m3) — r° = jr? + m3c2y/r2 + m3c2. Square again: 
4 
c 
z (mt — mz — m3)? — ne (my —mi- m3) +74 = ES (m3 + m3) + mmác* 
4 
C 
P C — m3 — m3)? — Amm] = r°? m$ + m$ + m? — m3 — m3) = r*mic* 


2 





c 
= z [mi + ms + m3 — 2m?m3 — 2m?m3 + 2mžm3 — 4m3m3| 
1 
c 
r= mí + m5 + m$ — 2m?m2 — 2m?m2 — 2m2m2 |. v 
1 
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Two-particle Decays 


pau? f lu (r)|28 (mc — u) Ż du 
— 8rňm J, vm u 


ma +m3)c 





Last integral sends u to m,c and hence r to: 
C 
To = m Ví + ms + mi — 2mm} — 2mim3 — 2mm 


ro is particular value of |p2| that is consistent with conservation of 
energy. In comprehensive notation: 


SIpl 
= ES VAL 
8r himic 


where |p| is magnitude of outgoing momentum, given in terms of three 
masses (eqn of ro). 
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Two-particle Decays 


For massive particles we have the very useful expression: 


SIp| 





IMI? 


= 8rhim?c 


Example decay: p > 1+1 





For massless particles using conservation of energy we have 
|p| = (m1c)/2. Thus decay rate becomes: 





[= IM]? 


1 a 


Example decay: n — y +y 
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Three-particle Decays 


@ In the previous two examples, we have shown how the phase 
space for the decay rate of a 2-body decay can be integrated 
completely without any information about M. 

e For 3-body decays (and beyond), this is no longer possible, as the 
amplitude will typically depend non-trivially upon several of the 
phase space integration variables, so that we have to do the 
integration by hand for each specific M. 
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Golden Rule for Scattering 


Golden Rule for Scattering 
Suppose particles 1 and 2 collide: 


14233+4+4+4---4+n 


The cross section is given by 


sh? 


-4 Sprm} — (mme) 


4 š 
Teri — m; fe) 8 (e) oi 


f L4] 2x)*8* (m + pa — p — Pn) 





Where the variables are all defined as in the decay case. 
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Golden Rule for Scattering 


The phase space is essentially the same as the decay case: integrate 
over all outgoing momenta, subject to the three kinematical constraints 
(every outgoing particle is on its mass shell, every outgoing 
energy is positive, and energy and momentum are conserved), 
which are enforced by the delta and theta functions. 

We can simplify cross section by performing the p? integrals: 


sh? 
o= ———— 
4V (pi ` p2)? — (m mc?)? 
= 1 dďp; 
(2m)? 


Í I-A A (2x)*8*(p + pa — ps: — Pn) 





x ——— 
F 2 2-2 
j=32 pj + mic 


with 


f= dil ende 
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2-body Scattering Example 


Consider the process 





Uz 


in the CM frame. If the amplitude is M, calculate the differential cross 
section. 


Py P2 6 
— © © e--L---- 


Before "d After 
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2-body Scattering Example 


and you can show (using E? = (pc)? + (mc?)? a few times) 


(p1 - P2)? — (mı mac?? = (Ey + Ez)Ipil/c 


Which allows us to write the cross section 


2 
Si?c fet Sintp-bn-r» ps — pa) d 


3 
2 Xe NY P3 d'p 
64r?(E1 + E2)lpıl p + mic /pž + mic? dd 


Delta function can be rewritten as: 


ETE 


&*(pi + pa — ps — pa) -sí A 


— p§ -) 5° (p3 + pa) 
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2-body Scattering Example 


Now the integral will pick out values corresponding to p4 = —ps. 


nM Sc f 2 
o =(—) = | |-41 
F3 (E + E)lpil 
5 C: + E)/c — ps + mie? — Jp + ie 
O o d 
V pi + mic? pi + mie 


Using spherical co-ordinates: 


d P; 


d?ps — r?drdQ 


where r = |p3| and dQ = sin8 d0 dý 
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2-body Scattering Example 


We can not carry angular integration as M? depends on direction of ps 
and magnitude. 


We can thus get differentical cross section as: 


2 oo 
do - (s) ao. | WAR 
dQ 81 (Ei + E2)|p; | 0 
5 | + E)/c— Jr + me — Jr? + "ie 
2 
X =... fy 
r? + mig |r? + mgo? 
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2-body Scattering Example 


The integral then yields (see the variable substitution as in case of 
two-body decay) 


Cross Section for 2-body Scattering 


do (ej SMĚ pi 
dQ 8n) (E + Eo)? |pil 
where |p;| is the magnitude of outgoing momentum and |p| is 
magnitude of incoming momentum. 





The dimensions of M are 
(mc)? 


where n is the number of external lines in the Feynman diagram. 
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The Feynman Rules for a Toy Theory 





OK, now we know how to take a known amplitude (M) and use it to 
calculate a decay rate or a cross section....but how do we get M?? We 
will start with Griffith's Feynman rules for a toy theory. Let's say we live 
in a world with only 3 kinds of particles (A, B, C) and only one vertex: 


B 


C 


This vertex is also the diagram for A — B + C decay (assume A is 


heavier than B and C combined). These three particles are all spin-O 
and each is its own antiparticle. 
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The Feynman Rules for a Toy Theory 


We want to be able to calculate the lifetime of A to lowest order and 
calculate the scattering processes A + A S B+B 





and A--B—A-B 
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The Feynman Rules for a Toy Theory 


@ Draw the Feynman diagram(s) for the lowest-order process, ie, 
with minimum number of vertices. There may be more than one. 
© Label the incoming and outgoing 4-momenta p4, po, ..., Pn. Label 


the internal momenta as q1, q»,..., qn. Put arrows on each line to 
indicate the "positive" direction. 


@ For each vertex write down a factor of 
-ig 
where g is called the coupling constant. For this toy g has 


dimensions of momentum, normally it is dimensionless. 
@ For each internal line, write a factor 


i 
2?  m?c? 
qt- mic 
where qj is the 4-momentum of the line and mj is the mass of the 
particle the line describes. 
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The Feynman Rules for a Toy Theory 


@ For each vertex, write a delta function of the form 
(27)*6* (Kk, + ko + kg) 


where the "k"s are the four-momenta coming into the vertex. This 
imposes conservation of energy and momentum at each vertex. 
Q For each internal line, write down a factor 


1 
— —d^gi 
(27)“ ai 
and integrate over internal momenta 
@ You will be able to write the result containing a delta function 


(27)*6*(pi + po +: — Pn) 


erase this factor. Multiply by i in order to get M. 
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The Feynman Rules for a Toy Theory 





To sum up the Feynman rules: 
The amplitude is given by: 
M 


i (vertex factors) 
(propagators) 
(momentum conservation) 
(integration over internal momenta) 
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Lifetime of the A 


The simplest possible diagram is that for the decay of A 





B 


C 


There are no internal lines, only one vertex so we get 
© vertex coupling: —ig 
© conservation of E and p: (27t)*5*(p, — pg — pc) 
@ erase the delta function. Multiply by i in order to get M. 


M-g 
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Lifetime of the A 


© From Fermi's Golden Rule, the decay rate is given by: 


_ _ Sip! IMP? 
8rhm?c 


in the rest frame of A. 
e With S = 1, m4 = ma, M = g and |p| representing the magnitude 
of the spatial momentum of either B or C, we find that 


_ _ ip! 
8rhm*,c 


and the lifetime is M~! 
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Scattering in a Toy Model 


Consider the lowest order scattering diagram in our toy model 
ue b 
B 





This has two vertices (each —ig) and one internal line: 
RUF M 
g = mée? 
The vertices give 2 delta functions 
(21)*8 (p: — ps — q), (27)*6^ (pa + q- pa) 
and the internal line gives an integration 
1 
(2n)* 





d*q 


Jyothsna (CHEP) Nuclear & Particle Physics Course October 2018 Lectures 349 / 431 


Scattering in a Toy Model 


We then have 
1 
-i(21)*g* f 50" (D: — Ps - q)ó* (po + q - pa)d“g 
q? - mec 


The second delta picks out the values at q = p4 — p» leaving 


np, 1 


-i 21)^6^(pi + p2 — ps — 
9 Ga poj m za Qr) (pi + pa — Ps — pa) 





Remove the remaining delta and we have 


g? 
© (Pa— p)? - mc? 
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Scattering in a Toy Model 


But wait! What about: 





That has the same initial and final state but a different assignment of 
final momenta (p3 © p4). The total amplitude (at lowest order) for this 
process is therefore 


2 2 


g g 
DTE a GM 
(pa — p2)? - mec? — (ps- p2)? - mec? 


We add the amplitudes of all diagrams with same same initial and final 
state before plugging into the cross section expression. 
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Higher Order Diagrams 





9 Of course, we have only looked at the lowest order ("tree level") 
scattering diagrams. What about diagrams with more vertices? 


e These diagrams will generate extra contributions to the amplitude: 


Ma B«c 





MA+B>C+D 


= gA,+9°A3+g°As +- 
= gb» g* Aa + gP7ls - - 


@ |f g << 1 then successive terms in the perturbation series become 
smaller and smaller corrections to the amplitude. 
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Corrections to A Decay 





There is one third-order diagram to consider: 


This O(g?) contribution to the O(g) tree-level diagram will come in to 
the decay rate in a coherent sum: 


IMÊ = IJA + As? 
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Corrections to A Decay 





We can also see diagrams for A decay (assuming A is sufficiently 
heavy) which do not result in the same final state: 


B B 
B C 


A 


This is A > 3B + C. The total decay rate for A could be something 
like: 


r(A > anything) = Ta sc +lassesc + l'ASBOCC +: 
2 
C; Du Masc] + C2 Y Ma->BBBC 


2 
+03 X Ma>eccc| 


where the C's are from the Golden Rule. This expression involves both 
coherent and incoherent sums. 
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Í 


Scattering Corrections 


We also have many possible diagrams to correct scattering. For 
example: 








five 'self energy’ (left), two ‘vertex correction’ (center), one 'box' 
diagram (right) 

seven more with external B lines twisted. 

In total there are 15 of these fourth-order diagrams. Clearly this could 
be a lot of work! 
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Scattering Corrections 


Let's just look at the “vacuum polarization" correction to the amplitude 








A RR 
See pages 218-219 of your text for details, but we end up with 


g i 1 [ d*q 
2n/ [(p:— ps)? - mz c?P J [(p1 — ps - q)? - má c?)(q? — mac?) 


Attempting to calculate this integral will give you fits. At large q we have 
oo 1 " 
f GO = Ing“ = c 
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M = i( 


Regularization 


9 The first step in dealing with loop integral divergences is called 
regularization. This is an artificial adjustment to the integral so 
that it can be solved (“sweeping the infinities under the rug"). 


e Introduce a cutoff mass. Introduce a factor: 
—M2c2 
(q? — M?c?) 
9 This cutoff mass M is taken to be some very large number. 
Effectively you change the upper limit of integration from co to M. 


© Some of the infinities introduced in the coherent sum over 
amplitudes now cancel....but not all. We now have a finite part and 
an infinite part. 


9 Now we need renormalization to rescue us. 
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Renormalization 


© All divergences in the final physical observables (o, I) appear as 
additions to the masses and coupling constants. 

9 So, we can define renormalized masses or couplings which 
absorb the divergences. We then assume it is the renormalized 
quantities that we have been observing all along. This means that 
the bare quantities aren't observable. 


Mprysical = M+6M 
Ophysical = 9+ ôg 

e We use the physical masses and couplings in the Feynman rules. 
(So, | guess we better measure them) 

e Once infinities are resolved we still have finite contributions from 
higher-order diagrams. These are functions of the 4-momentum of 
the line in which the loop is inserted. So, effective masses and 
coupling constants depend on energy - they run! 

© Renormalizability is a feature of all QFTs. 
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Renormalization 











JFCARTIER. 


ALRIGHT RUTH, I ABUT GOT THIS ONE RENORMALIZED, 
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